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ABSTRACT 

The  objective  of  this  analytical  study  was  to  develop  a  rapid 
theoretical  analysis  on  approximate  half  wavelength  elements  of  an 
antenna  configured  in  a  co-planar,  symmetrical  array.  A  computer  pro- 
gram was  written  employing  the  method  of  moments  approach  to  the 
solution  of  Hall  en's  integral  equation  with  an  approximate  two  term 
entire  domain  expansion  assumed  for  the  current  distributions.  With 
the  solution  of  the  current  distributions  on  each  element,  additional 
calculations  were  made  for  the  input  impedance  and  admittance  values, 
field  distributions,  power  gain,  and  a  graphical  output  of  the 
radiation  pattern. 
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I.   INTRODUCTION 

A.  DEVELOPMENT  OF  THE  STUDY 

A  goal  in  the  theory  of  multi -element  antenna  design  is  to  have  a 
tool  which  easily  analyzes  the  behavior  of  an  array.  The  equations 
that  describe  the  behavior  of  a  one  voltage-driven  element  with  a 
simple  sinusoidal  current  distribution  are  readily  available  in  many 
texts  [1].  However,  upon  the  introduction  of  a  second  or  third  element 
with  different  physical  dimensions—either  voltage  driven  or  parasitic-- 
the  behavior  without  the  assumption  of  a  simple  sinusoidal  current 
distribution  becomes  rather  complex  and  lends  itself  to  a  computer 
oriented  solution. 

A  large  portion  of  the  theory  employed  in  the  programming  of  this 
software  package  was  presented  and  discussed  in  a  seminar  at  the 
University  of  Mississippi,  sponsored  by  the  National  Science  Foundation 
in  1972  [2]. 

The  computer  program  of  this  study  was  written  in  the  XDS  Fortran 
IV  language.  It  was  specifically  to  be  used  in  the  computer  laboratory 
of  the  electircal  engineering  department  at  the  Naval  Postgraduate 
School,  Monterey,  California  [3]. 

B.  ANALYSIS  OF  THE  PROBLEM 

A  brief  synopsis  of  the  theory  of  this  study  is  as  follows: 

1 .  The  Development  of  Hall  en's  Integral  Equation 

2.  The  Solving  of  Hall  en's  Equation  for  the  Current  Distribution 
by  the  Use  of  an  Assumed  Two  Term  Approximate  Solution  Employing  the 
Method  of  Moments 


3.  The  Development  of  the  Associated  Equations  for  the  Relative 
Vertical  Plane  Pattern,  Relative  Horizontal  Plane  Pattern,  and  the  Gain 
of  the  Array 

4.  An  Outlining  of  the  Computer  Programming  Procedures 

C.  OUTLINE  OF  THE  RESULTS 

An  antenna  array  of  from  one  to  ten  elements  was  modeled  in  free 
space.  The  simple  physical  parameters  of  the  designed  array  are  the 
input  to  the  computer  program.  All  calculations  were  made  by  the 
computer,  producing  an  output  of:  (a)  a  listing  of  the  current 
distribution  on  each  element  of  the  array;  (b)  the  input  impedance 
and  admittance  quantities  for  each  of  the  voltage-driven  elements; 

(c)  a  listing  of  the  relative  field  distributions  about  the  array; 

(d)  a  graphics  output  of  the  radiation  patterns;  and  (e)  the  values 
of  the  gain  in  magnitude  and  in  decibels. 

D.  SCOPE  AND  LIMITATIONS 

The  procedures  outlined  were  completed  for  a  symmetrical,  co-planar 
array  of  from  one  to  ten  elements  in  free  space  with  each  element  being 
approximately  a  half  wavelength  long.  No  boom  effects  were  included. 
Each  of  the  ten  elements  may  be  voltage-driven  or  an  array  of  one 
voltage-driven  element  and  nine  parasitic  (reflectors  or  directors) 
elements  may  be  configured.  The  latter  includes  possible  Yagi-Uda 
arrays.  The  elements  may  be  spaced  arbitrarily  apart  and  the  center- 
driven  voltage  may  be  any  complex  value.  The  radii  of  the  elements  are 
restricted  to  the  thin  wire  approximation  (a  <  .01  a).  The  graphical 
radiation  pattern  outputs  were  of  the  relative  vertical  and  horizontal 
planes. 


No  optimal ity  procedures  were  introduced.  Further  study  in  this 
area  could  be  extended  by  introducing  a  finitely  conducting  ground 
plane  to  the  model;  include  elements  which  are  other  than  co-planar; 
or  include  elements  which  are  not  symmetrical  about  the  axis,  but  are 
at  angles  to  the  axis. 
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II.  FORMULATION  OF  HALLEN'S  INTEGRAL  EQUATION 

A.  THE  DIFFERENTIAL  EQUATION  FOR  THE  VECTOR  POTENTIAL  FROM  THE 
ELECTRIC  FIELD  INTENSITY 

The  basis  of  this  study  is  the  calculation  of  the  approximated 
two-term  current  in  the  solution  of  Hallen's  integral  equation,  as  all 
other  quantities  are  determined  from  this  data.  In  this  section,  a 
brief  discussion  of  the  formulation  of  Hallen's  integral  equation  is 
presented.  The  objective  is  to  obtain  an  expression  of  the  current 
distribution  of  a  cylindrical  center-fed  voltage  driven  antenna  in 
terms  of  its  length  and  diameter. 

Consider  a  voltage  driven  element  as  depicted  in  Figure  1.  It  is 
known  that: 

e;   =   ez  (n-i) 

E'   =  E  (II-2) 

p     P 

In  this   case,   the  electrical    field  intensity,  F,  and  the  vector 
potential,  A,  exist  only  in  the  Z-component.     Therefore,   the  equation, 

.   2 

E     =  ^—  v(v-A)    -  ju,A  (II-3) 

becomes , 

— J-  +  k%  =  jk2r7  (n-4) 

dr  L  — 

CO 

Equation   (II -4)   is  correct  everywhere,  where  JMx,y,z)   is  the  vector 
potential   and   E"z(x,y,z)   is  the  electric  field  intensity.     This   equation 
is  solved  on  each  element  of  an  array  of  co-planar  elements   in  the  YZ 
plane  as  seen  in  Figure  2.     The  radii   of  the  elements   in  this  figure 
are  represented  as    'a   '. 
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FIGURE  1.  COMPONENTS  OF  THE  ELECTRIC  FIELD  AT  THE 
SURFACE  OF  THE  ANTENNA 
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Evaluating  (II-4)  on  the  surface  of  each  of  the  N  elements 
produces : 

-   y    +  k2¥w>  ■    z  n  n 

dZ  (o 

Ze(-hn,  +hn)  (II-5) 

Solving  (II-5)  with  Bn  and  Cn  being  arbitrary  constants  and  V   being  an 
arbitrary  point  on  the  'nth'  element  gives: 

Az  (Vyn>z)  =  Cncos(kZ)  +  Bnsin(kZ) 

+^-  f         E7(a  ,y  ,S)  sin(Z-S)k  dS 
5  n 

n  -  1,  ....N  (11-6) 

simplifying, 

Av(a   ,y   ,z)     =     C  cos(kZ) 
Z v   n ' J  n  *   '  n       v      ' 


^/c-n  EZ(WS>   sin(Z"S>k 


S=0 

n     =     1,    ....N  ( 1 1  - 7 ) 

Equation  (1 1-7)  is  achieved  with  the  requirements  that  the  array  of 
elements  are  perfectly  symmetric  in  the  YZ  plane  so  that  the  scalar 
potential  is  zero  at  the  center  of  each  element. 

Now  an  expression  for  the  vector  potential,  A,  will  be  expressed 
in  terms  of  the  antenna  current  and  the  two  expressions  will  be  set 
equal  to  each  other,  for  the  formulation  of  Hall  en's  integral  equation. 
The  only  unknowns  of  the  equation  will  be  the  current  distribution  and 
a  constant  of  integration. 
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B.  THE  DIFFERENTIAL  EQUATION  FOR  THE  VECTOR  POTENTIAL  FROM  ELEMENT 
CURRENTS 

At  an  arbitrary  point,  'P',  in  Figure  2,  the  vector  potential  can 

be  written: 

N 
A(?)  =  E  M?)  (II-8) 

i=l  ] 

where  A^(r)  is  the  vector  potential  due  to  the  current  on  the  ' i th ' 

element  and  is  expressed  as: 

h 


A.(F)  =  /—  f     i  I.(zl) 


e-jk|r  -  r! 


dz! 


(H-9) 


zr~hi 


r  -  r. 


where, 


r  -  v- 


x2  +  (y  -  yi)2  +  (z  -  z!)2 


(11-10) 


In  the  above  equation,  x,  y,  and  z  are  the  coordinates  of  the  point 
P'.  The  y'  and  z.'  are  the  coordinates  on  the  'ith'  element. 
Substuting  ( 1 1 -9 )  and  (11-10)  into  ( 1 1-8)  produces  (11-11). 


A  (?) 

zv  ' 


h  i—   —  i 

N  r   i        -jk  r  -  r'-\ 


1=1   1   1 


r  -  r. 


Equation  (11-11)  is  the  vector  potential  at  an  arbitrary  point  in 

space  due  to  al 1  N  antenna  currents.  Evaluating  (11-11)  along  the 

surface  of  the  'nth'  element  produces  (11-12). 

h 

(H-12) 


where 


A2(an,yn,z) 


Rri  " 


N  -h. 


-hl        R'ni  dzi 


a2  +  (y  -  y!)2  +  (z  -  z!)2 
n    VJn   Jv     v     i' 


(H-13) 
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C.  HALLEN'S  INTEGRAL  EQUATION 

Setting  equation  ( 1 1-7)  and  equation  (11-12)  equal  to  each  other 
along  each  of  the  N  elements  results  in  Hallen's  integral  equation 
(11-14). 


.-JkR^(z.Z') 


4lT  i=i*/z'=-hi  n     Rni(z,z') 


dZ 


Cncos(kZ)  +  ^ f      Ez(an,yn,S)  sin(Z-S)k  dS 

n  =  1 ,  . . . .  N 


(11-14) 


where 


Rni(Z,Z')  - 


a2  +  (y„  -  yO2  +  (Z  -  Z')2 


(H-15) 

Since  the  current  in  each  element  is  an  even  function  of  Z,  the 
limits  of  the  integrals  can  be  changed  to  zero  and  h-.  With  this 


change  the  set  of  N  equations  in  (11-14)  becomes  (11-16),  where, 

-JV, 


yz) 


60 


^  sin(kZ) 


t(z) 


— -  cos(kZ)  =  ~  cos(kZ) 

JkRni(Z,Z')   e-JkRni(Z,Z') 

+  -R-xzrz-r 


Kni(Z'r)  =  "  KrfU.ll 


(H-17) 
(11-18) 
(II-19) 


With  this  formulation  of  Hallen's  integral   equation,  it  was  this 
study's  objective  to  solve  the  N  set  of  equations   (11-16)   for  the 
unknowns  of  I-.    ....    I*.  and  the  constants  of  integration  C-j    ....   C^,  as 
a  function  of  the  arbitrary  distance  along  each  element  Z'. 
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III.  THE  TH EORY  OF  THE  TWO  TERM  APPROXIMATE  SOLUTION  TO 
HALLEN'S  INTEGRAL  EQUATION 

A.  TWO  TERM  APPROXIMATION 

Even  though  the  N  set  of  equations  (II--16)  may  look  rather  compli- 
cated, there  exists  only  two  sets  of  unknown  quantities.  Knowing  the 
values  of  the  parameters  of  an  array  of  elements,  the  currents  and 
constants  of  integration  are  easily  determined. 

The  technique  of  the  two-term  solution  assumes  that  the  current  on 
the  element  may  be  expressed  as  (III— 1). 


I1(Z)  =  E  bj  sin  pir/2h.  (h.  -  Z) 


p=l 


Expanding  (III-l)  produces  (II 1-2) . 


(III-l) 


I/(Z)  =  b]  sin 


ir/2hi  (hi  -  Z) 


b?  sin 


2ir/hi  (hi  -  Z) 


(IH-2) 


In  equation  (III-l),  the  'i'  represents  the  number  of  the  element 
and  the  ' p'  is  the  number  of  the  coefficient  in  the  sinusoidal  expan- 
sion. This  expansion  sets  the  number  of  unknowns  at  two  current 
coefficients  and  one  constant  of  integration  per  element.  Therefore, 
there  are  three  unknowns  per  element  in  the  set  of  equations  (11-16). 
Applying  (III -2)  to  (11-16)  produces  the  set  of  N  equations  in  (II 1-3) . 
For  an  example,  a  three  element  array  (N=3)  would  have  nine  unknowns 
in  the  set  of  three  equations.  They  are,  as  seen  in  (III-3),  b-j  ,  b^, 
b-,  ,  b?,  b,  ,  bp,  C-, ,  C2,  and  C^. 

B.  METHOD  OF  MOMENTS 

In  order  to  determine  these  unkowns,  the  method  of  moments  is 
employed  [4],  The  set  of  testing  functions  selected  are  a  set  of  delta 
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functions,  which  enables  the  usage  of  the  point  matching  procedure.  In 
this  study,  the  selected  number  of  matched  points  per  element  equal  the 
number  of  unknowns  per  element,  that  is,  three.  The  point  matching 
process  then  produces  three  equations  for  every   one  of  the  N  equations 
in  (II 1-3) .  Each  integral  in  (III -3)  must  be  solved  for  Z  equalling 
one  of  the  designated  matched  points  and  V    existing  as  the  dummy 
variable  of  integration.  The  notation  used  in  (II 1-3)  is  as  folic 
h. 


Smp  "/  1  sin  ^2hi  <hi  "Z') 


K.ni(Znni5Z')  dZ'       (III-4) 


Z'=0 

Here,  the  * i '  represents  the  number  of  the  element,  ' p '  is  the 
number  of  the  coefficient  in  the  sinusoidal  expansion  (in  this  study 
p=l  or  p=2),  and  'm'  is  the  number  of  the  matched  point  on  a  particular 
element,  'n'.  Employing  this  notation  and  arranging  all  the  unknowns 
to  the  left  hand  side  in  (II 1-3 ) ,  produces  (II 1-5),  where  there  are 
three  equations  for  e\/ery   one  in  (III-3).  In  (III-5)  as  previously, 
the  unknown  quantities  are  the  b  coefficients  and  the  C  constants  of 
integration.  Consider  again,  the  example  of  the  three  element  array 
(N=3)  with  its  nine  unknowns.  Equation  ( 1 1 1 -5 )  becomes  nine  simultane- 
ous integral  equations  with  nine  unknowns.  Every  S   becomes  a  solvable 
integral  with  the  selection  of  the  matched  points  (Z  ). 

C.  MATRIX  FORMATION 

Rearranging  (II 1-5)  into  matrix  form  produces  (II 1-6 ) ,  which  is  in 
the  form  AX=B.  The  "A"  matrix  must  be  filled  in  with  the  desired 
number  of  zeroes,  depending  on  the  number  of  elements,  N.  The  'A' 
matrix  will  always  be  square  and  be  a  3N  by  3N  size  matrix.  The  'X' 
and  'B'  matrices  will  always  be  a  3N  by  1  size  matrix.  Clearly,  the 
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'A1  and  'B'  matrices  contain  all  known  quantities  and  the  'X'  matrix 
contains  all  the  unknown  quantities. 
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IV.  ASSOCIATED  EQUATIONS  OF  THE  STUDY 

A.  ADMITTANCE  AND  IMPEDANCE 

The  solving  of  the  current  distribution  from  the  matrix  equation 
(III -6 )  is  the  essence  of  this  study.  From  the  current  distributions, 
the  admittance  and  impedance  quantities  are  found  by  using  the  terminal 
current  (current  at  center  of  element,  where  the  voltage  is  applied) 
and  Ohm's, Law. 

B.  RELATIVE  VERTICAL  AND  HORIZONTAL  PLANE  PATTERNS 

The  calculations  begin  with  the  solving  of  the  differential 
equation  (IV-1)  in  the  far  field  region  [5]. 

dH   =  3(1)  =in  9  dZ  (UM) 

9         cS  A 

The  representation  of  this  equation  is  found  in  Figure  (3).  The 
relative  vertical  plane  pattern  EQ  in  free  space  is  found  from  (IV-1) 


by  equation  (IV-2) . 


EQ  =  ZH   -  120ttHa  (IV-2) 

e     9       9 


The  current  distribution  calculated  previously  is  pictured  in 

Figure  (4).  The  approach  used  in  solving  (IV-1)  was  by  superposition 

of  the  two  currents.  The  current  i]  may  be  substituted  into  (IV-1)  as: 

.  ju(t  -  s/c) 

i]  sin  [k(L/2  ±  Z)]  e  (IV-3) 

and  the  current  I„  as : 

jw(t  -  s/c) 
\\  sin[2k(L/2  i  Z)]   e  (IV-4) 
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Point 


FIGURE  3.   RELATIONS  FOR  SYMMETRICAL,  THIN,  LINEAR,  CENTER-FED 

ANTENNA 
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FIGURE  4.  TWO  TERM  EXPANSION  OF  THE  CURRENT 
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In  Figure  (3),  the  distance  V  to  a  distant  point  is  found  to  be: 

S  =  r  -  Zcose  (iv-5) 

Substituting  (IV-3)  and  (IV-5)  into  (IV-1)  produces  an  integral 
equation  of  the  form  in  (IV-6). 


/ 


ax  . 


ax 


e  sin(c  +  bx)dx  =  -~ «  [a  sin(c  +  bx)  -  b  cos(c  +  bx)]  (IV-6) 

a^  +  bd 


Carrying  through  the  integration  and  making  use  of  the  trigonometric 
identity  .in  (IV-7)  and  the  relationship  that  h  =  L/2  from  Figure  (3), 


the  result  is : 
H 


cf>l  2iTr 


?  2 

sin  e=l-  cos  ( 


cos(khcose)   -  cos(kh) 
sine 


(IV-7) 


(IV-8) 


where  I*'-  ije^  "  r/c) 

Relating  this  to  (IV-2)  and  setting  rrl  gives: 


Ee1   =  J60I] 


cos(khcose)   -  cos(kh) 
sine 


(IV-9) 


Employing  (IV-4)  and  (IV-6)  in  (JV-1)  produces  the  second  components 
of  the  superposition. 


H 


jig 

<t»2  :  27F 


2sine 


4  -  cos' 


cos(khcose)  -  cos(kh) 


sine 


(IV-10) 


Relating  (IV-10)  to  (IV-2)  and  setting  r=l  produces  (IV-11). 
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J60I 


r 


2sine 


4  -  cos  e 


cos(khcose)  -  cos(kh) 


sine 


(IV-11) 


Now  for  multiple  symmetrical  elements,  a  phase  term  is  introduced  due 
to  the  separation  of  elements  along  the  Y  axis.  In  taking  this  into 


27 


account,  the  time  term  of  I       is  dropped  and  with  the  substitution  of 

■; " 
k=u/c  and  r  =  yisino,the  form  of  I       becomes  that  of  (IV-12). 

, ,         ,     jky.sine 
Ip     =  Ip  e       1  (IV-12) 

Finally,  substituting  (IV-12)    into  (IV-11)  and   (IV-9),  the  far 

field  pattern  in  the  Y-Z  axis   is  achieved,   using  the  superposition 

equation   (IV-13). 

En     =  Efll   +  En0  (IV-13) 

e         el        e2  ' 

The  relative  horizontal   plane  pattern  as  depicted  in  Figure  (5)   is 

found  with  equation   (IV-14). 

,         ,       jky.costj) 
E^     =  j60(l]   +  V?)   e       1  (IV-14) 

C.      POWER  GAIN 

The  power  gain  was  calculated  for  comparison  with  an  isotropic 
source  [6].     The  gain  is  equal    to  the  ratio  of  the  power  intensity  to 
the  power  density.     The  gain  is  then  expressed  as   in   (IV-15). 

G  =  ML  (IV-15) 


w 
in 


where  from  the  Poynting  vector, 


120* 

With, 


r2  ? 

W   =  t4vt  |E|  (IV-16) 


Win=    |I|2Rin  (IV-17) 


The  gain  with  r~l   is  expressed  as   in  (IV-18) 


in2 
G  =  Hi (IV-18) 

30|I|2R.n 
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FIGURE  5.  CONFIGURATION  OF  PARALLEL  ELEMENTS 
IN  THE  XY  PLANE 
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V.  PROGRAMMING  PROCEDURES 

A.  ANALYSIS 

With  the  equations  of  the  unknown  available,  the  writing  of  the 
software  which  would  produce  computer  solutions  with  graphical  outputs 
began.  The  program  consists  of  the  main  program,  one  function  sub- 
program, and  four  subroutines.  The  following  is  a  brief  outline  of 
the  procedures  used,  in  the  writing  of  the  program. 

B.  OUTLINE  OF  STEPS  EMPLOYED 

1.  The  number  of  arrays  to  be  analyzed,  1  to  99,  and  the  designated 
number  of  the  graphical  output  device  [7]  to  be  used  is  read  in. 

2.  The  number  of  elements  in  the  array  to  be  analyzed,  1  to  10, 
and  the  wavelength  in  meters  is  read  in. 

3.  The  applied  complex  voltage,  value  of  'Y1  along  the  Y-axis  in 
meters,  half  length  of  element  ' h1  in  meters,  and  radius  'a'  in  meters 
is  read  in  for  each  element. 

4.  Constants  are  formed  including  the  designation  of  the  position 

of  the  matched  points  for  the  solution  of  the  integral,  S  . 
v  s     nip 

5.  That  portion  of  the  'A'  matrix  involving  the  integral  S   is 
solved.  The  integration  routine  employs  Weddle's  Rule  [8]  which 
approximates  the  function  with  a  sixth  degree  polynomial. 

6.  The  remainder  of  the  'A'  matrix  is  then  completed  by 
appropriately  filling  in  the  zeroes  and  the  function,  t(Z  ). 

7.  The  'B'  matrix  is  completed  by  solving  the  function,  V  (Z  ). 

8.  With  the  'A'  and  'B'  matrices  formed,  the  'X'  matrix  is  solved 
for  by  using  a  Gaussian  Elimination  routine  [9]  including  the 
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decomposing  of  the  'A'  matrix,  the  solving  for  the  'X'  matrix,  and 
then  the  attempt  to  improve  the  'X'  matrix  for  the  most  accurate 
solution  possible. 

9.  With  the  solution  of  the  'X'  matrix  achieved,  the  values  of 
the  current  are  outputed  starting  at  the  center  (Z=0)  and  for  every 
tenth  of  the  distance  along  each  element  until  the  end  of  the  element 

10.  For  every   element  which  is  voltage-driven,  the  input  impedance 
and  admittance  values  are  outputed. 

11.  Starting  at  theta  equal  to  zero  degrees,  and  for  every   degree 
through  360  degrees,  the  values  of  E  -,  and  E  ~  for  each  element  are 
summed  together. 

12.  These  values  are  normalized  for  output  and  the  appropriate  XY 
values  are  calculated  for  the  graphical  computer. 

13.  Starting  at  phi  equal  to  zero  degrees  and  for  every   degree  to 
360  degrees,  the  values  of  E.  for  each  element  are  summed  together. 

14.  These  values  are  normalized  for  output  and  the  appropriate  XY 
values  are  calculated  for  the  graphical  computer. 

15.  The  XY  coordinates  of  the  relative  patterns  to  be  drawn  are 
sent  to  the  graphics  computer  for  display. 

16.  The  power  gain  in  magnitude  and  decibels  is  computed  and 
outputed. 
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VI.  CONCLUSIONS 

A.  COMPARISON  OF  RESULTS 

The  parameters  of  two  simple  antenna  configurations  were  used  as 
input  data.  Results  of  this  study,  were  then  compared  for  degrees  of 
accuracy. 

B.  EXAMPLE  ONE.  TWO  SYMMETRICALLY  VOLTAGE  DRIVEN  ELEMENTS 

The  first  antenna  consisted  of  two  identical  symmetrically  voltage 
driven  elements  with  a  wavelength  set  equal  to  one  meter  (300MHz). 
Each  element  was  driven  with  the  complex  voltage  of  1+jO.  The  separa- 
tion between  the  two  elements  was  .01  meters  (y, =-0.005,  y?=0.005). 
Both  elements  were  0.50  meters  in  length  (L=2h).  The  radius  of  each 
element  was  0.000277  meters.  Computer  results  produced  by  the  program 
of  this  study  for  the  current  distribution  are  displayed  in  Figure  6. 
These  results  show  good  correlation  with  those  of  King's  first-order 
theory  [10]  and  a  two  term  sinusoidal  approximation  as  depicted  by 
Butler  [11]. 

With  this  study's  calculated  current  distributions,  the  input 
admittances,  input  impedances,  relative  vertical  pattern,  relative 
horizontal  pattern  and  power  gain  in  magnitude  and  decibels  were  also 
calculated  as  shown  in  Table  I. 

C.  EXAMPLE  TWO.  SIMPLE  SYMMETRIC  CENTER  FED  DIP0LE 

The  second  example  is  a  simple  symmetric  center-fed  dipole,  with  a 
wavelength  set  equal  to  one  meter  (300MHz).  The  complex  voltage 
applied  was  1+jQ.  There  is  no  separation  involved  in  this  example 
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IMAGINARY  PART 


REAL  PART 


-7.5  -5.0 


KING'S  FIRST-ORDER  THEORY 


•o  o  e  •  •   BUTLER  TWO-TERM  SINUSOIDAL  BASIS 
+  +  +  ■+   RESULTS  OF  THIS  STUDY 


FIGURE  6.  CURRENT  DISTRIBUTION  FOR  TWO  SYMMETRICALLY 
VOLTAGE  DRIVEN  ELEMENTS 
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REAL  CURRENT 

IMAGINARY  CURRENT 

DISTANCE  ALONG  ELEMENT 

(ma.) 

(ma.) 

FROM  CENTER  (meters) 

5.28236974 

-2.37995399 

.000 

5.16780772 

-2.47847437 

.025 

4.92962568 

-2.50660181 

.050 

4.57696180 

-2.45519573 

.075 

4.12109759 

-2.31881755 

.100 

3.57492582 

-2.09652088 

.125 

2.95246973 

-1.79220618 

.150 

2.26848157 

-1.41511770 

.175 

1.53813553 

-0.97857729 

.200 

0.77681740 

-0.50012839 

.225 

0.00000000 

-0.00000000 

.250 

INPUT  ADMITTANCE  REAL 

INPUT 

ADMITTANCE  IMAGINARY 

(mi li Siemens) 

i 

(mi li Siemens) 

5.28236974 

-2.37995399 

INPUT  IMPEDANCE  REAL 

INPUT 

IMPEDANCE  IMAGINARY 

(ohms) 

(ohms) 

157.36507827509 

70.90030883439 

OMEGA          RELATIVE  VERT.  PATTERN 

RELATIVE  HORZ.  PATTERN 

(cleg.) 

(absolute) 

(absolute) 

000 

1.00000000000 

.99950656036 

010 

.97786160215 

.99952143819 

020 

.91417041942 

.99956427764 

030 

.81632464159 

.99962991265 

040 

.69429231999 

.99971022788 

050 

.55865218871 

.99979611278 

060 

.41741211316 

.99987663247 

070 

.27633644475 

.99994227433 

080 

.13729063227 

.99998611978 

090 

.00000000000 

1.00000000000 

100 

.13729063227 

.99998511978 

110 

.27633644475 

.99994227433 

120 

.41751211316 

.99987663247 

130 

.55865218871 

.99979611278 

140 

.69429231999 

.99971022788 

150 

.81643564159 

.99962991256 

160 

.91417041942 

.99956427764 

170 

.97786160215 

.99952143819 

180 

1.00000000000 

.99950656036 

The  gain  of  the  array  is  3.09  =  4.90db. 

TABLE  I.  COMPUTER  OUTPUT  FOR  TWO  SYMMETRICALLY  VOLTAGE 

DRIVEN  ELEMENTS 
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(y-|=0.0).  Following  the  length  and  radii  restrictions,  the  length  of 
the  element  was  0.?5  meters  and  the  radius  was  0.0069  meters. 

Computer  results  produced  by  the  program  of  this  study,  are  dis- 
played in  Figure  7.  These  results  are  plotted  along  with  a  two  term 
sinusoidal  approximation,  and  a  five  term  sinusoidal  approximation 
as  displayed  by  Butler  [12]  and  measured  data  by  Mack  [13].  The 
additional  data  calculated  as  listed  in  Example  one  is  shown  in  Table  II 

D.  RECOMMENDATIONS 

More  complicated  arrays  were  also  analyzed  by  the  computer  program 
of  this  study.  These  outputs  yielded  results  which  were  closely 
related  to  those  of  other  two  term  approximations  for  current 
distributions  [14]. 

In  Examples  one  and  two,  the  matched  points  used  were  Z=0.0, 
Z=0.5h-,  and  Z=h..  It  is  recommended  in  a  study  by  Darko  Kajfez  [15] 
that  for  the  closest  comparable  results,  using  the  two  term  approxima- 
tion of  this  study,  the  matched  points  should  be  located  at  Z=0.0, 
Z=0.75h.,  and  Z=h. . 

This  study  found  that  the  number  of  increments  that  the  integral 
function  is  divided  into  when  being  approximated  by  a  sixth  degree 
polynomial  employing  Weddle's  Rule,  will  in  some  cases  have  an 
unwieldy  effect  in  the  solutions  of  the  integrals.  It  is  recommended 
that  the  integer,  75,  be  employed  in  the  determination  of  the  number  of 
increments. 

Recommendations  for  further  study  are  to  employ  a  different  complex 
numerical  integration  routine  with  the  theory  of  this  study.  Possibly 
additional  programming  could  be  written  in  order  to  allow  interactive 
graphics  for  design  and  testing  procedures. 
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RESULTS  OF  THIS  STUDY 


FIGURE  7.  CURRENT  DISTRIBUTION  FOR  SIMPLE  SYMMETRIC 
CENTER-FED  DIPOLE 
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MAGNITUDE 
(ma) 

10.06945176 
10.31137742 
10.29630974 
9.99144960 
9.37425727 
8.43725169 
7.19133077 
5.66716441 
3.91449351 
1.99942153 
0.00000000 


ARGUMENT 

DISTANCE 

ALONG  ELL" 

(deg.) 

FROM 

CENTER 

(meters) 

20 

.91494045546 

.000 

23 

.88716574162 

.025 

26 

.55377940962 

.050 

28, 

.92611960589 

.075 

30 

.97158254089 

.100 

32 

.68483783763 

.125 

34, 

,07240413432 

.150 

35, 

.13865883043 

.175 

35, 

,89288284490 

.200 

36, 

,34223309019 

.225 

36, 

,49146393058 

.250 

INPUT  ADMITTANCE  REAL 
(mi li Siemens) 


IMPUT  ADMITTANCE  IMAGINARY 
(mi li Siemens) 


9 

.40598988 

-3.59460891 

INPUT  IMPEDANCE  REAL            INPUT 

IMPEDANCE  IMAGINARY 

(ohms) 

(ohms) 

92 

.76685977913 

35.45193909900 

OMEGA 

RELATIVE  VERT.  PATTERN 

RELATIVE  HORZ.  PATT 

(deg.) 

(absolute) 

(absolute) 

000 

1.00000000000 

1.00000000000 

010 

.97740371048 

1.00000000000 

020 

.91249376631 

1.00000000000 

030 

.81307695801 

1.00000000000 

040 

.68974122014 

1.00000000000 

050 

.55322994009 

1.00000000000 

060 

.41223256076 

1.00000000000 

070 

.27215639036 

1.00000000000 

080 

.13498559880 

1.00000000000 

090 

.00000000000 

1.00000000000 

100 

.13498559880 

1.00000000000 

110 

.27215639036 

1.00000000000 

120 

.41223256076 

1.00000000000 

130 

.55322994009 

1.00000000000 

140 

.68974122014 

1.00000000000 

150 

.81307695801 

1.00003000000 

160 

.91249376631 

1.00000000000 

170 

.97740371048 

1.00000000000 

180 

1.00000000000 

1.00000000000 

The  gain  of  the  array  is  1.64  ::  2.15db 
TABLE  II 


COMPUTER  OUTPUT  FOR  SIMPLE  SYMMETRIC 
CENTER-FED  DIPOLE 
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APPENDIX  A 
GUIDE  TO  OPERATION  OF  PROGRAM 

Below  is  a  step  by  step  procedure  to  follow  in  order  to  operate  the 
computer  program,  at  the  Naval  Postgraduate  School  computer  laboratory. 
The  software  requires  the  use  of  the  digital  computer  (XDS9300)  and 
the  graphics  computer  (AGT). 

A.  FORMATION  OF  THE  DECK 

1.  The  first  card  is  a  'BOOT'  card  which  is  found  in  the  card 
stalls  in  the  laboratory. 

2.  The  second  card  is  the  'JOB'  card  which  is  followed  by  the 
cards  of  the  program  through  to  the  'DATA'  card  as  shown  in  the 
enclosed  computer  program  printout. 

3.  Following  the  'DATA'  card  is  the  initial  data  card  which 
contains  three  digits  in  the  first  three  columns  of  the  card.  The 
first  two  digits  denote  the  number  of  arrays  to  be  analyzed  (01  to  99). 
The  third  digit  denotes  the  assigned  number  of  the  graphics  computer 

to  be  used  (1  or  2) . 

4.  After  the  initial  data  card,  the  input  data  of  each  of  the 
arrays  is  formed.  Each  array  will  have  one  card  plus  as  many  cards  as 
there  are  elements  in  the  array.  The  first  card  of  each  array  will 
have  the  number  of  elements  in  the  array  in  the  first  two  columns 

(01  to  10).  In  the  next  twenty  columns,  the  value  of  the  wavelength 
in  meters  is  to  be  entered. 
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5.  Following  the  'wagelength'  card,  there  is  one  card  for  each  of 

the  elements.  The  format  for  each  of  these  cards  is  as  folio*  . 

Columns  01  to  15:  Real  value  of  complex  voltage  applied  to  element. 

Example  1.0 

Columns  16  to  30:  Imaginary  value  of  complex  voltage  applied  to 

element.  Example  -1.0 

Columns  31  to  45:  Value  of  the  Y-axis,  indicating  the  position  of 

the  element  in  meters.  Example  0.0 

Columns  46  to  60:  Value  of  half  length  of  element  (h)  in  meters. 

Example  0.25 

Columns  61  to  75:  Value  of  radius  of  element  in  meters. 

Example  0.001 

6.  As  an  example,  in  order  to  run  example  one  and  two  of  section 
VI  of  this  study,  the  cards  after  the  'DATA'  card  would  appear  as 

f ol 1 ows  : 

021 

022 

02   1.0 

1.0  0.0  -.005  0.25  .000277 

1.0  0.0   .005  0.25  .000277 

01   1.0 

1.0  0.0  0.000  0.25  .0069 

B.  PREPARATION  OF  XDS9300  DIGITAL  COMPUTER 

1.  Press  IDLE  switch  at  the  XDS9300  control  panel. 

2.  Hold  the  RESET  switch  depressed  and  press  the  POWER  switch  in 
order  to  energize  the  computer. 

3.  Turn  switch  at  teletype  to  the  'on'  position. 

4.  Press  the  READY  button  on  the  line  printer. 


Card 

1: 

or 

Card 

1: 

Card 

2: 

Card 

3: 

Card 

4: 

Card 

5: 

Card 

6: 
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C.  PREPARATION  OF  THE  AGT  GRAPHICS  COMPUTER 

1.  Press  ON  at  the  operator's  control  panel  (OCP). 

2.  Wait  5  seconds,  then  press  RESET  on  the  OCP. 

3.  Press  the  ON/START  switch  at  the  disk  drive  of  the  AGT  to  be 
used. 

4.  Press  HALT,  RESET,  and  BTSP,  in  that  order  on  the  OCP. 

5.  Verify  that  the  'location  counter'  (behind  the  door  of  the 
left-most  rack)  indicates  zero  or  all  lights  out.  If  not,  press 
RESET  and  BTSP  again. 

6.  Load  a  paper  tape  marked  DISK  BOOT  LOADER  (located  behind  the 
door  of  the  left-most  rack)  into  the  tape  reader  of  the  AGT's  teletype, 

7.  Verify  that  the  knob  on  the  front  of  the  AGT's  teletype  is  in 
the  'line'  position. 

8.  Put  the  tape  reader  switch  to  'start.' 

9.  Now  verify  that  the  'location  counter'  indicates  110  (octal). 
If  it  is  not,  return  to  step  1  of  this  section  and  begin  the 
procedure  again. 

10.  When  step  9  is  satisfied,  verify  that  the  switches  in  the 
center  of  the  OCP  have  the  configuration  24000  77776  (octal). 

11.  Press  HALT,  RESET,  RUN,  and  PULSE  1  on  the  OCP. 

12.  If  the  AGT's  teletype  does  not  respond  now  with  a  request 

for  the  operator  to  input  MO/DA/YR,  press  HALT,  RESET,  RUN,  BUTTON  29, 
and  PULSE  1.  If  teletype  does  not  respond  now,  return  to  step  1  of 
this  section. 

13.  Type  RESET( "GATED" ,101 ) !  on  the  AGT's  teletype. 

14.  Type  GATED!  on  the  AGT's  teletype. 

15.  Turn  vertical  and  horizontal  gain  knobs  fully  clockwise. 
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D.  READING  IN  PROGRAM  DECK 

1.  Place  cards  in  card  hopper. 

2.  Push  POWER  ON  button,  at  the  card  reader. 

3.  Push  START  button,  at  the  card  reader. 

4.  Push  IDLE,  RESET,  RUN,  and  CARDS  on  the  XDS9300  control  panel 
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